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Solutions for linear conservation laws with gradient constraints 


Jose Francisco Rodrigues Lisa Santos 


Abstract 

We consider variational inequality solutions with prescribed gradient constraints for first order lin¬ 
ear boundary value problems. For operators with coefficients only in L 2 , we show the existence and 
uniqueness of the solution by using a combination of parabolic regularization with a penalization in the 
nonlinear diffusion coefficient. We also prove the continuous dependence of the solution with respect 
to the data, as well as, in a coercive case, the asymptotic stabilization as time t —> +oo towards the 
stationary solution. In a particular situation, motivated by the transported sandpile problem, we give 
sufficient conditions for the equivalence of the first order problem with gradient constraint with a two 
obstacles problem, the obstacles being the signed distances to the boundary. This equivalence, in special 
conditions, illustrates also the possible stabilization of the solution in finite time. 


Dedicado a Joao Paulo Dias, no seu ativo septuagesimo aniversario! 


1 Introduction 

Several works have developed solutions u = u(x , t) to the linear equation of first order 

d t u + b ■ Vu + cu = /, (1) 

for t > 0 and x in an open subset f l of 'R N , where b = b(x,t) is a given vector field and c = c(x,t ) and 
/ = f{x,t) are given functions. 

The well-known DiPerna and Lions theory of renormalized solutions, when b is given in Sobolev spaces, 
has been extended by Ambrosio to BV coefficients for the Cauchy problem and has found several applica¬ 
tions in the study of hyperbolic systems of multidimensional conservation laws (see, for instance I|, for an 
introduction and references). The initial-boundary value problem for (|T|) with a C 1 vector field b has been 
studied in the pioneer work of Bardos jj] using essentially a L 2 approach for the transport operator. This 
method also holds for Lipschitz vector fields, as observed in [8], and was extended by Boyer [5j for solenoidal 
vector fields in Sobolev spaces that do not need to be tangential to the boundary of Cl, i.e. b ■ n ^ 0 on dCl 
for t > 0. 

The delicate point is then to prescribe the boundary data to the normal trace of b on the portion of the 
space-time boundary T_ C dCl x (0, T) where the characteristics are entering the domain Qt = Cl x (0,T). 
In the case when T_ does not vary with t , Besson and Pousin [3J have treated the initial-inflow problems for 
the continuity equation ([I]) with L°° velocity fields b with c = V-6 = divfo also in Recently Crippa 

et al. E have also considered this problem without that restriction on T_ and with similar assumptions on 
b in BV. 

Here we are interested in the initial-boundary value problem for 0 under the additional gradient con¬ 
straint 

\S7u(x,t)\<g{x,t), (x, t) e Q t , (2) 

where g = g{x,t ) is a given strictly positive and bounded function. This problem was already considered in 
m in the framework of a quasilinear continuity equation 

d t u + X7 ■ $(u) =F(u) (3) 
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and a Lipschitz semilinear lower order term F = F(x,t,u), with a gradient bound in (J2]) that may depend 
also on the solution but not on time. As observed in ED], in the linear transport equation 0, corresponding 
to 

<&(u) = bu and F(u) = f + (V ■ b — c)u 

with regular coefficients and g = g{x) independent of f, the problem is well-posed in terms of a first order 
variational inequality with the convex set 

I 9 = {u £ Hq(£1) : |Vu(x)| < g(x) a.e. x £ O}. (4) 

In |20j it is also proved the existence and asymptotic behaviour of quasivariational solutions for positive 
nonlinear gradient constraints g = g(x,u) depending continuously on the solution u = u(x,t). Here H^fY) 
denotes the usual Sobolev space of functions vanishing on the boundary dfl, as the gradient bound allows 
to prescribe values on the whole boundary. Moreover, it allows also to consider the data 6, c and / only in 
L 2 (Qt), provided c — | V b is bounded from below. 

A motivation for the constraint (|2|) applied to the equation (|T|) is the “transported sandpile” problem. 
Following Prigozhin p3J[15], the gradient of the shape of a growing pile of grains z = z(x,t) characterized 
by its angle of repose a > 0 is constrained by its surface slope, i.e. g = arctana. A general conservation of 
mass, in the form § with $ = —/uS/u + bu and source density F, with transport directed by b and dropping 
flow directed to the steepest descent — /zVu, should be then subjected to the unilateral conditions 

g, > 0, |Vu| < g and |Vu| < g => ji = 0. 

We illustrate this problem with the interesting example of the one dimensional special case announced in 
{19] : = (0,1), b = 1 = g , i.e. a = \ and f(x,t) = t. Taking as initial condition the parabola zq(x) = — |® 2 , 

up to the point £o = V3 — 1, and the straight line zq(x) = x— 1, for £o < x < 1, the profile of the “transported 
sandpile” growth attains a steady state exactly at t = This happens with the first free boundary point 
£(£) increasing from (o up to f = ! 2 , touching then the boundary x = 1, and decreasing till the midpoint 
x = |. At this point, the free boundary £(f) meets a second increasing free boundary £(£) = 2(f — 1), that 
appears at t = 1 and increases up to the final stabilization at t = |. 



Figure 1: The free boundary of the transported sandpile problem at t = 0,3/4,9/8 and 5/4. 
The explicit sandpile profile is given by 


tx — 1 X 2 

if 0 < x < £(£) and 0 < t < 1, 

x — 1 

if £(£) < x < 1 and 0 < t < \. 

1 — X 

if £(f) < x < 1 and | < t < 1, 

X 

if 0 < x < C(f) and 1 <t< |, 

tx — \x 2 

if £(f) < x < £(t) and 1 < t < 

x — 1 

if £(t) < x < 1 and 1 < t < |, 

\ - k- \\ 

if t > 


z(x,t ) = < 
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where £(f) = t— 1 + \J{\ — t ) 2 + 2, if 0 < £ < and £(t) = t + 1 — \/(i + l) 2 — 2, if | < t < 

It is clear that z(t) € Kj C 

We introduce the function d(x) = \ — \x—\ \ and the convex set 

Ky = {d£ Hq(LI) : —d{x) < v(x ) < d{x) a.e. x £ (0,1)} D Ki. 

Since d t z + d x z = t in A = {{x,t) £ Qt ■ \d x z(x,t)\ < l}, by simple computation and integration in Qt , 
we easily conclude that z, which (using V = sup and A = inf ) can be written as 

z(x, t) = ( — d(x)) V ((tx — l;x 2 ) A d(x )), 

is then the unique solution in Ky of the variational inequality 


/ (d t u + d x u — t)(w — u) > 0 V w{t ) £ K(), 0 <t <T, u(0) = z 0 . 

J q t 


( 5 ) 


But since zo,z(t) £ Ki, z is also the solution of the variational inequality ([5]) with w(t) £ Kj C Ky, which 
has at most one solution also in the convex set Ki, defined as in 0 with < 7 = 1 . 

In Section 2 we establish the existence and the uniqueness of the solution of the first order variational 
inequality associated with the general linear equation |l]) in a family of time dependent convex sets with 
gradient constraints of the type Q with g = g(x,t). We improve the results of (20] under general square 
integrability assumptions on the coefficients and on the data, by direct estimates in the parabolic-penalized 
problem and passage to the limit, first in the penalization parameter £, and afterwards in the regularization 
parameter 6. The continuous dependence of the solution with respect to the gradient constraint variations 
in L°° , to the coefficients of the operator and the data in L 1 , is proven in Section 3 under the weak coercive 
condition (J7| , as well as the asymptotic convergence towards the unique stationary solution under the stronger 
coercive assumption (231. 


Finaly, in Section 4, we consider the special case of a constant vector b, with <7 = 1 and / = f(t) bounded, 
to show the equivalence of the variational inequalities with the gradient constraint and with the two obstacles, 
i.e. with the signed distances to the boundary constraints on the solution. This is a first result of this type 
for first order variational inequalities, similar to the elliptic well-known case of the elastoplastic torsion 
problem (see, for instance, [TO] and its references) and to the parabolic case without convection considered in 
Earns, where it was shown that this equivalence is not always possible in the general case. With additional 
conditions, that include the above one dimensional transported sand pile problem, we establish the finite 
time stabilization of the solution. This extends to the convective problem a similar result by Cannarsa et 
al. Bj and raises the interesting open question of establishing more general conditions on the finite time 
stabilization of evolutionary problems with gradient constraints. 


2 Existence and uniqueness of the variational solution 

Let fl be a bounded open subset of with a Lipschitz boundary dfl and, for any T > 0, denote Qt = 
fl x (0, T). 

Assume that 

b£L 2 (Q T ) and cGL 2 (Q t ), (6) 

and there exists l £ R such that 

c - • b > l in Q t , (7) 

being this inequality satisfied in the distributional sense, since V • b does not need to be a function. 

In addition we also suppose given 


( 8 ) 


with 


/ £ L 2 {Qt) and uq £ K g (o)) 
g £ W 1 ’ oo (0,T;L oo (n)), g > m > 0. 


( 9 ) 





4 


Solutions for linear conservation laws with gradient constraint 


As in Q, we define, for t > 0, 

K 9 ( t ) = {y € ifg(f2) : |Vv(ir)| < g(x,t) for a.e. in x G U}. 

Consider the following variational inequality problem: To find u, in an appropriate space, such that 
' u(t) € K fl ( t ) for a.e. t £ (0, T), u( 0) = uq, 


[ d t u(t)(v - u(t)) + I b(t) ■ Vu{t){v - u(t)) + f c(t) u(t)(v - u(t )) > [ f(t)(v - (10) 

Jq Jq Jq Jq 

Vv £ K g ( t ), for a.e. t £ (0,T). 

Theorem 2.1 With the assumptions problem ( |10[ ) has a unique solution 

u £ T°°(0,T; W 0 1 ’ oo (fl)) n^(Q T ), d t u £ L 2 {Q T ). 

Proof To prove the uniqueness of the solution we assume there exist two solutions u± and U 2 ■ Using 
U 2 = U 2 {t) as test function in (10) for the variational inequality of iq and reciprocally, setting u = iq — U 2 
at a.e. t > 0, we obtain 


j d t u(t)u(t) + j b(t) ■ Vu(t) u(t) + j c{t)u 2 {t) < 0. 
J Q J Q J Q 

Using |7]), for any v € 'if 0 oc (O), we have 

^ f b(t) ■ Vu 2 + f c{t)v 2 > l f v 2 
2 J n J n J n 


and, by approximation in of u(t), we obtain, 


j j \u(t)\ 2 +2l I \u(t)\ 2 < 0. 


By Gronwall’s inequality, we conclude u = 0 from u(0) = 0. 

To prove the existence of a solution, we consider a family of approximating quasilinear parabolic problems 
for u es , with e,S £ (0,1), defined as follows 


d t u eS - <5V • (fc £ (|Vu E,5 | 2 - g 2 )'\7u eS ) + b 5 • Vu El5 + c 5 z t E<5 = f 5 in Q T , 
« Ei =0onafix(0,T), (11) 

m ei5 (0) = Mq in Uq, 


where b s , c* 5 , f s and ng are appropriate regularizations of 6, c, / and uo, respectively, with |Vitg| < g(0) 
and k e is a smooth real function such that k e (s) = 1 if s < 0 and k E (s) = e = if s > e. Notice that this 
problem has a unique solution u eS £ H 1 (0, T; L 2 (U)) fl-L°° (0, T; ilg (Q)) n^(Q T ), by the classical theory of 
parabolic quasilinear problems (see, for instance, HD]). 

We prove first several a priori estimates. 


Estimate 1 

\\k e {\Vu s5 \ 2 -g 2 )\\ L , {QT) <- 5 C x , 

for some constant C x dependent only on to, \l\, ||/||l 2 (q t ), |MU 2 (q t ) and IKIUqn)- 

Multiplying the equation of the problem ( |TT| ) by u eS and integrating over Q t = Ux]0,t[, we have 

\ f \u £ \t)\ 2 + S [ fc £ (|Vu E<5 | 2 — < 7 2 )| Vw E<5 | 2 + [ (b 5 ■ Vn £i ) u eS + [ c VT = [ fu sS + J [ 
1 Jn JQ t J Q t J Qt J Qt 1 Jn 


( 12 ) 


,e |2 
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Observing that 

[ (b 5 ■ Vu sS ) u e5 = [ (y ■ b 5 ) (u e5 ) 2 , 

JQt 2 JQt 

and using the coercive inequality for the regularized coefficients 


c s - ■ b s = (c - ■ b) * ps > l * ps = l, 


we have 


[ \u^(t)\ 2 + 5 [ k e (\X7u eS \ 2 - g 2 )\S7u eS \ 2 < \\f S \\ LH Q T) \\u^\\ LHQt) + h\u 0 \\ 2 L2{n) + \l\\\u^\\ 2 L2{Qt) . 

J £1 J Q+ ^ 


Hence, by the integral Gronwall’s inequality, there exists a positive constant Ct, independent of e and 
5, such that 

IK 5 !U 2 (Qt) < c t 

and so 

Sf fc e (|Vu ei5 | 2 — g 2 )\Vu e5 \ 2 < C', 

JQt 

where C' = C'(||/|U 2 (Q r ), ||u 0 || L 2 (n)) , \l\). 

On the other hand, we observe that 

[ k e (\\/u eS \ 2 — g 2 )\X/u eS \ 2 = f k e (yu sS \ 2 -g 2 )(yu s5 \ 2 -g 2 )+ [ k e (\Vu eS \ 2 - g 2 )g 2 . (13) 

J Qrp J Qrjp J Qrjp 

Since k e (s) = 1 for s < 0 and k E {s)s > 0, for all s > 0, then 


[ k e (\S/u e5 \ 2 — g 2 ){\S/u eS \ 2 — g 2 ) = f 


{|Vu e 5 | 2 <(f 2 } 

r 

{\Vu^\ 2 >g 2 } 


fce(|Vu e5 | 2 — 5 ' 2 )(|Vw e<5 | 2 — g 2 ) 

k s (yu eS \ 2 - g 2 )(\Vu eS \ 2 -g 2 )>-[ g 2 . (14) 

JQt 


From (13) and (141 we obtain 


(15) 


M|VH 2 - g 2 ) < ^ (^J k e {\X7w\ 2 -g 2 )\X7w\ 2 + g 2 ^ < Qc" + llslli^Q,,)) < ^C x , 

where C x = Ci(m, ||/|| l 2 ( q t ), \\gh 2 (Q T ), |KIU 2 (n), \l\). 

Estimate 2 

l|V« e5 || inQT ) < D S , 

where, for any <5 > 0 and any 1 < p < oo, the constant D$ depends only on p, m, |/|, ||/||l 2 (q t )> ||'?I|l 2 (q t )) 
||wo||_L2(n) and a negative power of 5. 

From © we know that 

[ k e (\Vu eS \ 2 - g 2 ) < \c u 

JQt S 

where C\ is the positive constant of Estimate 1. So, 




M |v« rf | 2 - s 2 ) = 


{\Vu eS \ 2 >g 2 +e} 


| Vu £ °\ z -g z 


'{| V u eS | 2 > g 2 +e} 


and recalling that, for all s > 0 and all j e N, e s > 4y, we get, for any j € N, 


f (|Vt* rf | 2 -g 2 Y<j\e> f 

J{\Vu‘ s \ 2 >g 2 +e} Jf 


iVu eS l 2_ 2 l 

e -«- < -JWC 1 . 


{|Vu ei 5 | 2 >g 2 +E} 
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Given 1 < p < oo, we have 

f \Vu eS \ 2p = 
J Q T 


l{\Vu £S \ 2 <g 2 +e} 


| W 


e8\2p 


' {|Vn e 5 | 2 >p 2 +e:} 


| W 


eS\2p 


(16) 


and, since g is bounded, we can estimate, for any p £ N, the second integral in the second term of (16) as 
follows, 


/ {\\'u eS \ 2 >g 2 +e} 


\S7u s5 \ 2p < 


l{\Vu^\ 2 >g 2 +e} 


3=0 W 


The first integral in the second term of (16) is clearly bounded since 

\Vu e5 \ 2p < [ ( g 2 + lY 

J Qt’ 


J {\ Vu £S \ 2 <g 2 +e} 

and the conclusion follows easily, first for 2 p € N and afterwards for any 1 < p < oo. 

Estimate 3 

ll^t M ''‘ 5 ||| 2 (Q T ) < 4(ll^lli a (Q T ) + ^3 He 5 II z , 3 (Q^)) ll^ 7M ‘' <5 |li«(Q T ) + (1^) 

where, for 2 < s < , C 3 is an upper bound, independent of q = s 2s 2 , of the Poincare constant for W 0 1,? (fl) 

and C 4 is a positive constant depending only on to, |Z|, ||/||l 2 (Q t ), llffll vvi,~( 0 ,T ; L“(n)) and ll u o|U 2 (n)- 

We multiply the equation of problem (11) by d t u eS and we integrate over Q t , noting that d t u eS = 0 on 

dfl x (0,T). Denoting </> £ (s) = / k e (r)dT , we have 

Jo 

[ \d t u eS \ 2 + S -[ ^AM\^ s \ 2 -g 2 ))+s f k e (\Vu sS \ 2 -g 2 )gd t g 

JQ t 1 J Qt at J Qt 


+ [ (b s -S/uAdtu 65 + f c s u e 5 d t u eS = [ f s d t u eS 
J Qt J Qt ^Qt 


We choose 2 < s < 2 _„ and q = _ s .,, and so we have - + - + \ = 1. Then 


I (b S -VuAd t uA < \\b 5 \\ L s iQT) \\Wu e 5 \\ Lq{ Q T) \\d t u £ 5 \\ L 2 {QT) < ||blL(Q T )l!V^|lL ( Q T) + ^||a t ^||i 2(0T) , 
J Qt 


Qt 
and 


So 


c 5 u eS d t u eS 


< l|c d |U.(Q T )||« ea |U<(Q T )||ft« ea |U* WT ) < V\\l*{Q t )\W \\l<,(Q t ) + -j\\ dtU 


j Q \d t uA 2 <\\f 5 \\ 2 L*(Q T ) + yje(\Vu s \0)\ 2 -g 2 (0)) S - Jj e (\Vu eS {t)\ 2 -g 2 {t)) 

+ Cl||3||L=o(Q T )||9 t fl , ||L“(Q T ) + (ll^lli^Qx) + ^J c 5 ||L(Q t ))II^ m ^IIl‘!(Q t )> 


being C q a Poincare constant. Observe that, since fl is bounded we may find a positive upper bound C 3 of 
C q , independently of q < 00 . 

On one hand 

[ 0 £ (|V^(O)| 2 -<? 2 (O)) <0, 
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because |Vm £i 5 (0)| = |Vit§| < g(0). On the other hand, if we set A = {{x,t) £ Qt ■ \Vu eS {x,t)\ < g{x,t)}, 
we have 


<f) e (\S7u sS {x, t )\ 2 — g 2 {x,t)) = \Wu eS (x,t )\ 2 — g 2 {x,t ) > —g 2 {x,t) for a.e. (x,t) £ A, 
4> e (\Vu eS (x,t )\ 2 — g 2 (x,t)) >0 >—g 2 (x,t) for a.e. {x,t) £ Qt\A. 

Consequently, for a.e. t £ (0, T), 

- / <t>e{\Vu eS {t )\ 2 ~ g 2 {t)) < ||sll!™ ( 0 iT;i 2 (n)) . 

J 17 

So, 

^l|dt ME< 5 || L 2 ( Q T) < ||/' 5 |||2 (Qr) + 2ll 5 , lll oo (o,T ; L 2 (n)) 

+ Ci|l5lU“(QT)ll^t5 , IU 0 °(i3T) + (II^IIl^Qt) + C?ll c<5 |lL a ((2T)) II^ 7 mE<5 |Il9(q t )5 

and the proof of Estimate 3 is concluded. 


By (15) and ( |T7| ), we know there exist constants D$, Cs and C 4 , independent of e, such that, for each 
N < p < 00, 

ll^ £ ^ll z>(o,t ; Wq , p (n)) - \\d t n sS \\ LHQT) < (||b 5 ||is(Q T ) + \\c 5 \\ls(Qt))Cs + C 4 . 

Since u eS is bounded in H 1 (0, T; L 2 (fl)) D L p (0, T; Wl' p {fi) fl independently of e £ (0,1), 

forp > N, by a known compactness theorem (|23j, page 84), {m £ i5 } £ is relatively compact in e >f([0,T]; e tf(fl)'). 
Then, at least for a subsequence, 

u eS —>■ u s in ffiQrp). 

e ->0 

The above estimates also imply that we may choose, always with fixed 5, 
u eS -^ u s weakly in L p (0,T-,Wq’ p LI)), 1 < p < 00, d t u eS -^ d t u s weakly in L 2 {Qt)- 


£—>-0 


£—>-0 


Given v £ L°° (0, T ; H g {fl)) such that v(t) £ K g ( t ) for a.e. t £ (0, T), we multiply the equation of problem 
(11) by v(t) — u e 5 (t), we use the monotonicity of k e and we integrate over fl x (s, t), 0 < s < t < T, to 


conclude that 

rt 


d t u e6 {v - u ed ) + 6 Vv VO - u ed ) 

Js JQ 

+ [ f b s ■ f7u eS {v-u eS ) + 
Js J 17 


r cu eS {v~u eS ) > f f f s (v-u eS ). 
17 J s j 17 


Letting e — > 0, since s and t are arbitrary, we obtain that for a.e. t £ (0, T) 


f d t u s (t){v(t) — u s (t)) + 5 f Vv{t) ■ V(u(t) — u s {t)) 

J n J n 

+ [ b\t)-f7u s {t)(;v(t) - u s {t)) + f c s {t)u s {t)(v{t) - u s {t)) > f f s {t)(v(t)-u s {t)), 

J 17 J 17 J 17 

for all v £ L°°(0, T; Hq {fl)), such that, v{t) £ K ff ( t ) for a.e. t £ (0,T). 

Set A e = {{x,t) £ Qt ■ \S7u eS {x, t)\ 2 — g 2 {x,t) > v 7 e}- Since fc e (|Vzt E,5 | 2 — g 2 ) > in A e , then we have 


14,1 = / 1 < / 

J A r J A 


k e {\S7u e °\ 2 — g 2 ) ^ C\ 

1 A r e ve ; 

ev-i 0 
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by (12), being C 1 a constant independent of e as we have seen. So we have 
[ (IW 5 ! 2 - g 2 ) + < liminf [ (\Vu sS \ 2 - g 2 - v / i) + 

Jq t £-i-0 Jqt 

= liminf / (\X7u eS \ 2 — g 2 — s/e) < lim Ms |AJ S = 0, (18) 

E ->- 0 JA e e ^° 

where Ms is an upper bound of || |Vu e5 | 2 — g 2 — sfe\\ l 2 (q t ), independent of e. Consequently, 

|V t/| < g a.e. in Qt 

and so u s (t) £ K g ( t ) for a.e. t £ (0, T). Let 2 £ L°°(0, T; Hq (11)) be such that z(t) £ K s (t)- Defining 
v = u + 9(z — u), 9 £ (0,1], then v(t) £ Using v(t) as test function in (10) and dividing both sides of 

the inequality by 9, we get 

f d t u 5 (t){z(t)-u s (t)) + S f Vu 5 (t)-V(z(t)-u 5 {t)) + 69 I \V{z(t) - u\t ))\ 2 
J fi Jd JCl 

f b\t ) • Vu 5 {t)(z(t) - u s (t)) + [ c s (t ) u s (t)(z(t) - u s (t )) > [ f s {t)(z(t) - u 5 (t)) 

Jq Jo, Jq 

and, letting 9 —> 0, we conclude that u s solves the following variational inequality 

u s (t ) £ K g ( t ) for a.e. t £ (0, T), u 5 (0) = «o> 


f d t u s {f)(v — u s (t)) + S ( Vu s (t) ■ V(v — u S (t)) 

J Q J Q 

+ [ b 5 (t) -Vu s (t)(v - u s (t)) + j c s (t)u s (t)(v-u 5 {t)) > f f s (t)(v - u s (t)) : 

JQ J Q 

\/v £ 3Cg( t ), for a.e. t £ ( 0 ,T). 


(19) 


Recalling the Estimate 3 we have 
\\9tn sS \\l HQT) < (||b' 5 ||i S (Q T ) ■+■ C 3 1 1 v || 


+ C 3 \\c°\\ 2 ,, (n ^)\\Vu sS " 2 


< {\\» 


<5 M 2 


LHQt) + C '3||C"||l,.(Q r )J 


51 | 2 


Qi 


(iv 


L-HQt) 
eS I 


u 


+ C 4 

- - v /e) + + / (5* + Ve 


' Qt 


Passing to the liminf when e —> 0 and arguing as in (18), we conclude that 


liminf [ (|Vu e5 |^ - g ^ - s/e) + = 0 

e ^° JQt 


and, consequently, 


\\dtU S \\ L 2(Q T ) < 4(||6 5 |||.(q t ) + C'3 ||c 5 |||»(q t ))||5|P j* . + ^4- 

L“ *(Qt) 


Observing that 

(II^IIl^Qt) + C 3\\c s \\ 2 L s { Q T ^)\\g \\ 2 2 _ 


L s 2 (Qt) s— ^2 


(l|bii 2( Q T) + c' 3 ||c 5 i| i2(Qr) )i| ff ||^ (Qr) , 


+ c 4 . 


we have the sequence {dtu s }s uniformly bounded in L 2 (Qt)- 

Moreover, the sequence is uniformly bounded in L°° (0, T; W 0 1,o ° (0)), independently of S, since 

each u s (t) belongs to So, there exists a function u £ L°° (0, T; W 0 1 ,oo (O)) D H 1 (0, T; L 2 (Mj) n c ^’(Q T ) 

and, at least for a subsequence, 


<5->0 
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-—- u weakly in L p (0, T; Wo’ p (Q)), 1 < p < oo, dtu s ——^ dtu weakly in L 2 (Qt )• 

< 5—^0 S —>-0 


Integrating in (19) between s and t, for 0 < s < t < T, and passing to the limit when <5 —> 0, we get 

n d t u[v — u) + f I fo • Vm(m — m) + f f cu(v — u) > f f f(v — u), 

for all v such that v[t) £ K g (t) for a.e. t £ (0 ,T). Since s and t are arbitrary, we can drop the integration in 
time. Since u 5 (t) £ K 9 ( t ) for a.e. t £ (0, T), the same holds for u(t), concluding that u solves the variational 
inequality ( 10 ). □ 


Remark 2.2 We observe that in the proof of the uniqueness of the solution it is sufficient to assume only 

b £ L 1 (Qt) and c £ L 1 (Qt'), 


instead of (| 6 |. 

Similarly, we may replace 0 by the different weak coercive assumption by assuming the existence of 
r£R, such that, in the sense of distributions, 


V • b > 


in Qt, 


in order to have also the uniqueness of the solution to the variational inequality (10). 

In fact, assuming that there are two solutions u\ and U 2 , we may choose for test function v = u± + 
( 2 S({u 2 — Mi) in the variational inequality for u\, where : R -> R is a sequence of C 1 increasing odd 
functions approximating pointwise the sign function sgn° and f is sufficient small. Then, choosing also 
v = U 2 + C 2 s<(mi — M 2 ) in the variational inequality for M 2 , we get 

/ d t u(t) s<;{u(t)) + / b(t) ■ Vu(t) s^(u(t)) + / c(t) u(t) s^{u(t)) < 0 
J Cl J Cl JCI 

Noting S^(s) = J s^(t) dr -—* |s| and rs^(r) —> |r|, by the dominated convergence theorem, we have 

4 [ K0I+ [ b(t)-V\u(t)\+ f c(i)|u(t)|<0 

al Jci JCl Jci 


and 


d 

dt 


\u{t)\+r [ |m(»| < 0. 

Jci 

Since m(0) = 0, by the GronwalTs inequality, we conclude the uniqueness from 


m\<e rt / |m(0)| = 0. 


3 Stability and asymptotic behaviour in time 


In this section, the stability of the solutions of the variational inequality (10), as well as its asymptotic limit 
when t —> +00 is based in the following Lemma, which is due essentially to [55]. 


Lemma 3.1 Fori = 1,2, let gi belong to L°°{Qt )■ If v 1 € L q (0 ,T;Wq’ p (CI)) , 1 < p, q < 00 , is such that 
Mi (t) £ K gi (t) for a.e. t £ (0 ,T) then there exists M 2 £ L 9 (0,T; hLg’^n)) such that M 2 (t) £ K g 2 (t) for a.e. 
t. £ (0, T) and a positive constant C such that 


II Ml v 2\\lv(0,T;Wo’ p (CI)) — C'lbl ff2||i“(Q T )- 
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Solutions for linear conservation laws with gradient constraint 


Proof Let a(t) = \\gi(t) - g 2 (f)||L°°(n)- Define ip(t) = 1 + C — and v 2 (t) = -^-rVi(t). 

v ' m yj[t) 

Since 

|V ® 2( ‘ )I = W« |V ”' ( ‘ )IS W) 3,{>} 


and 


9i(t) 


m 


:9i{t) < 92(t) 


ip(t) m + a(t ) 

then v 2 (t) £ for a.e. t £ (0, T). The conclusion follows immediately from 

1 


|V(vi - v 2 )| = 


1 - 


V’(t) 


|V«i| < ■ ^ \\gi - 92 \\l°°(Qi 


□ 


The continuous dependence result is a consequence of the boundedness of the solution and of its gradient, 
when we impose the weakly coercive assumption (J7|. 


Theorem 3.2 Fori = 1,2, let Ui denote the solution of the variational inequality (101 with data 
gi,uoi) satisfying assumptions ©-©• Then there exists a positive constant C = C(T), depending on T, 
such that 

||wi - fo2|||~(0,T;.L=(fr» < C(IKl “ u 02|| ^ 2 ( 0 ) + ||&l - MlfoQ-r) + 11 C 1 - C 2 ||li(Q t ) 

+ ||/l - /2 Hi , 1 (Qt) + llfll “ 52||i”(Q T ))- 

Proof Let u 2 be defined as in Lemma O for the solution u± and u\ be the corresponding function for u 2 . 


Using ii\ as test function in the variational inequality (10), we obtain 

[ dtui(i)(tii(t) - u 2 (t)) + f b x (f) • Vui(t)(Mi(t) - u 2 (t)) + [ ci(t)Mi(t)(u 1 (t) - u 2 {t)) 

J f2 J £2 J Q 

< f fi(t)(ui(t) - u 2 (tj) + [ ( d t ui{t ) + bi(t) • Vui(t) + Ci(t)ui{t) - /i(f)) (ui(t) - u 2 (t)) 

J Q JQ 

and a similar inequality is true using the variational inequality of u 2 , by replacing the data /i,i>i,Ci by 
f 2l b 2 , c 2 and u\ by u 2 . Then we have 

f d t (ui(f) - u 2 (t)) (ui(t) - u 2 (t )) + [ bi(t) ■ V(ui(t) - u 2 (t)) (ui(t) - u 2 {t)) 

Jn Jn 

+ [ Cl (t) (ui(t) - u 2 {t )) 2 < 0(f), (20) 

Jn 

with 

0(f) = / (d t ui(t) + bi(f) • Vui(t) + ci(t)ui(t) - fi(t)) (ui(t) - u 2 (t)) 

Jn 

+ [ (d t u 2 (t) + b(t) ■ Vu 2 (t) + c(t)u 2 (t) - f 2 {t)) (u 2 (t) - Ui(f)) 

Jn 

+ [ ( bi(t ) - b 2 (f)) • Vu 2 (f) (ui(f) - u 2 (f)) 

Jn 

+ [ [(ci(f) - C 2 (f)) u 2 (t) + (/i(f) - f 2 (t))] (tii(f) - w 2 (f))- 

Jn 

Using the boundedness of the solutions Ui , * = 1,2, and their gradients and recalling the L 2 (Qt) estimates 
of we have 

f 0(r)dr < C M (\\9i ~ 9 2 \\l°°(q t ) + ll b i “ MzfoQ-r) + II c i ~ c 2 ||li(q t ) + ||/i - MIl^Qt)), 

J o 
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where Cm is a positive constant depending on T, on the norms of the solutions and their derivatives (which 
can be bounded in terms of the data) and on the constant C of Lemma 3.1. 

Setting w = ui ~ U 2 in the inequality (20 1 , we obtain using Q , 

j t Jjw(t)\ 2 <2\l\ +2Q(t). 

Applying Gronwall’s inequality, we conclude 

I \Mt) - u 2 (t )\ 2 < e 2 |/|T (j|ui 0 - U 2 o\\ 2 L 2 {(1 ) +2G m (||5i ~ 92 \\l<*>(q t ) 

+ \\bl - b 2 \\L 1 (Q T ) + ||ci - C 2 \\ L 1(Q T ) + ll/l - /2 Hl 1 (Qt))^ * 


□ 


In order to consider the corresponding time independent solution to the first order variational inequality, 
we give stationary data foo , goo , boo , Coo satisfying the assumptions 


3oo G ffoo > m > 0, foo G Z^fi), 

booGL'ffl), CooGl 1 ^), 


Coo - • boo > A > 0 in 12 


in the distributional sense, where we set accordingly 

K goo = {w € Hq(LI) : |Vu>| < goo a.e. in 11}. 

Then, the stationary problem can be written as 

I boo * ^rioo (wJ Woo ) T / Coo Woo (w Woo ) I foo (w^ Woo) ; V W G 

J £1 J Q, J Q, 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 


Since the convex set K goo is bounded in Hq(Q) and the first order linear operator in the left hand side 
of (25) is pseudo-monotone, by the classical theory (see, for instance, [12]) it has a solution, which is unique 


by the strict coerciveness induced by the condition A > 0 in (231. 


In order to study the asymptotic convergence of the solution of the variational inequality (10 1 to the 


stationary solution of (25), we consider solutions global in time. This is easily obtained if we assume that 
are satisfied for any T > 0 and replace ©by 


g G W 1 ’ 00 (0, oo; L°°(f2)), g > m > 0. 


(26) 


We need an auxiliary lemma. 

Lemma 3.3 (0, pg. 286) Let (p : (0,oo) —> R. be a nonnegative function, absolutely continuous in any 
compact subinterval o/(0,oo), 4) £ Lj oc ( 0, oo) a nonnegative function and p a positive constant, such that, 

<p'{t) + pp>{t) < 4>(t), V t > 0. (27) 


Then, for any s,t > 0, 


ip(t + s) < e + 


1 

1 — e~v 





□ 
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Solutions for linear conservation laws with gradient constraint 


In order to apply this Lemma to 


tp(t) = ( | u(t) 
Jn 


- u„ 


t > 0 , 


we shall require the additional assumptions on the coefficients and on the data 

b £ L°°(0,oo; L 2 (Cl)) and c, / £ L°° (0, oo; L 2 (Cl)). 


(28) 


(29) 


Theorem 3.4 Assume that f,g,b,c,uo satisfy the assumptions (| 6 |-(| 8 |, ( [26] ) , ([29]) and /oo, < 7 oo, hoo, Coo sat¬ 
isfy the assumption (21) ,(22) and (23). Suppose, in addition, that 


pt -\-1 n pt -\-1 n rt -\-1 n 

/ / |/( T ) ~ foo\drdx —► 0, / / |6(r) - b^drdx —► 0, / / (c{t) - c^drdtx 

J t J J t J J t J 1*2 


and there exists 7 > \, such that, for some constant D > 0, 

hit) -ffoo|Uoc (n) < —, t> 0 . 


(30) 


If u and Moo are, respectively, the unique solutions of the variational inequalities (10) and (25) then, for 
every a, 0 < a < 1, 

u(t) —> itoo in ff 0 ' 01 (Cl) 

t—> OO 


Proof First we need to return to the estimate {l7| ) of the existence proof in order to prove that, under the 
additional assumptions of this theorem, there are positive constants A,B , independent of T, such that, 

||dt u lli 2 (nx(o,T)) < AVt + B. 

Since \X7u(x,t)\ < g(x,t) for a.e. (x,t) £ Q 00 = Cl x (0, 00 ) and g £ L°°(<3oo), we have now u £ 
L°°(0, 00 ; IF 1 ’ 00 (f2)). This yields the estimate 


I £ 2 (Qt) 


M < c g T 


where the constant c g > 0 is independent of T. Using similar estimates for ||/||^ 2 (q t ) with the constant 
c g replaced by c f = ||/||i 2(0iOO;L3(n)) , as well as for c b = || 6 |li 2 ( 0 iOO;I , 2 (n)) and c c = ||c||| 2(0oo;L2(n)) , we 
may conclude that the constant C\ = C\(T) of ( [12] ) , in the Estimate 1, grows also linearly with T, i.e. 
Ui < Co + c\T, where cq depends only on uq and c\ depends on to, c/, c g , c b and c c . Using this fact in the 
Estimate 3, we may now easily deduce ([3| from with s = 2 and q = 00 , since C 4 , depending on / and 
on Ci grows also linearly with T. 


Using Lemma 3.1 we choose u^ £ K 9 (t), for a.e. t £ (0,T), as test function in (10). Then 


in 


d t u(t)(u(t) - Uq 


/n 


b(t) ■ Vu(t)(u(t) — m 0 


In 


c(t)u(t)(u(t) - Moo) 


< [ f(t)(u(t) - Moo) + [ ( d t u(t ) + b(t) ■ Vu(t) + c(t)u(t) - f(t )) (uoo - Moo). 

dn Jn 


in Jn 

Analogously, with u(t) £ K Soo , for a.e. t £ (0,T), we obtain the inequality 


/ hoc 1 i^j(t) Mqo) “b j CoqU^q^u)!) Mqo) 

Jn Jn 

/ /oo (fj(t) Moo) “b / (hoo * ^ ^00 + Coo Mqo /oo) ( u(t ) u(t)) . 

Jn Jn 
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Then, simple algebraic manipulations lead to 


I df Zioo) (^(t) Woo) + / boQ ^oo) (w(t) Woo) T / Coo (w(t) Woo)(w(t) Wqo) 

Jo JO JO 

<0(t), (31) 

where 

0 (O= / (d t u(t) + b(t)-'Vu(t) + c(t)u(t)-f(t)) (Uoo-Uoo) + / (boo-Viroo+CocMoo-Zoo) (w(t)-M(t)) 

JO JO 

+ [ (b(t) - boo) • Vu(t)(Moo - w(t)) + [ (c(t) - C 00 )w(t )(« 00 - lt(t)) + [ (f(t) - foo)(u(t) - Moo), 
jf} 


Using (231 and the definition (28 1 , from ( |31| ), we obtain the differential inequality with p = 2A and where, 
taking into account Q, we may choose $(i) > 2 |©(t)| given by 

$(t) = C((H\/t + -B + C')Hg(t) - 3 oo||l~( 0) + II&W — ^oolliqo) + ll c (t) — Codl^qn) + ||/(t) — /oolUqo))- 


Then, using the assumptions and observing that the number 7 in (30) is greater than |, we have 


pt -\-1 rz-\-L 

/ $(r)dr <C (||/(r) - /oo||z,i ( o) + ||b(r) - boolkqo) + ||c(r) - || L i (n) )dr 

Jt Jt 

f t+1 1 

+ C' (t = +1)||3(t)- 3oo|U- ( o)rfT —> 0 . 

J ^ £—>-+oo 

w(t) 


rt+l 


Therefore, by Lemma 


3.3 


£—^-|-oo 


in L 2 (Q). Since u belongs to L°° (0, 00; TU 1 , 00 (fl)), the compact 


inclusion of W 1 , 00 (fl) in C' 0 ,Q: (Q) implies, first for a subsequence, and after for the whole sequence, that 

□ 


u(t) —> Uoo in C 0 ,a (f 2 ), concluding the proof. 
£—>•+00 


4 Finite time stabilization in a special case 

In this section we assume that dfl is of class 'if 2 and 

bet^, c = 0, <7 = 1 , z 0 eKi and f £ L°°(0,T). 

We consider the following two obstacles problem 

{ z(t ) £ Ky for a.e. t £ ( 0 ,T), z( 0 ) = w 0 , 

f d t z(t)(v — z(t)) + I b • V z(t)(v — z(t)) > f f(t)(v - z(t)), \/v £ Ky, for a.e. t £ (0, T), 
Jo Jo Jo 


(32) 


(33) 


where 

K() = {w £ Hq(Q.) : —d(x) < v(x) < d(x) for a.e. x £ U}. 

Here d(x) = d(x,dLi) is the distance function to the boundary dLl. Notice that d £ IUq , 00 (H), |Vd(x)| < 1, 
a.e. x £ Li and Ad < C for some constant C = C(fl) > 0. Observe that zq £ Ki C Ky. 


Theorem 4.1 Under the assumptions (32), the inequality ( |33| ) has a unique solution 

z £ L°°(0,T; Wq 1 ' 00 ^)) 0 tf 1 (0, T; L 2 (H)) ntf(Q T ), 


which satisfies |V.z| < 1 a.e. in Qt and is the unique solution of the variational inequality (10). 








14 


Solutions for linear conservation laws with gradient constraint 


Proof For e,5 £ (0,1), we consider the following family of penalized problems for z eS , 
( d t z eS - SAz eS + b ■ Wz eS + f (z eS - {z eS Ad) V {-d)) = f s in Q T , 
\ z eS ( 0) = z E 0 on n, z e 6 = 0 on dQ x (0,T), 


(34) 


where f s and z g are regularizations of the functions / and Zo, with |Vzq| < 1. This problem has a unique 
solution z eS £ Ff 2 , 1 (Q r ) ; since the operator 

(P £ v,w) = - [ (v — (v A d) V (— d))w (35) 

£ Jn 

is monotone (see, for instance, [T 2 ]h 

We obtain firstly an estimate of |V,s e ‘ 5 | on dfl x (0, T). Since dfl is of class if 2 , there exists r > 0 such 
that, if B r (x) denotes the ball with centre in x and radius r, then for all Xq £ dVL there exists y 0 £ such 
that B r (y 0 ) n = {a’o}. Placing the origin of the coordinates in the point yo, let i] £ (s) = e~ 7? and 

p(x) = d(x) + Me( 1 - y e (\x\ - r)), p(x) = -d(x) - Me( 1 - rj e {\x\ - r)), 


where M is a positive constant, depending on <5, to be chosen later. We show that p is a supersolution of 
(34). Analogously, it can be verified that p is a subsolution. We start by observing that 


p(xo) = 0 = z eS (xo,t) and p > 0 = z eS on dfl x (0, T). 


Since Zq £ K(), then 
We compute 


ip(x ) > d(x) > Zq(x). 


and 


Let 


d Xi ip(x) = d Xi d(x) + My/er} e (\x\ - r)^ and |Vy?| < 1 + My/e, 
dljp{x) = d^..d(x) — Mr) s (\x\ - r) ^ + Myfer} s (\x\ -r)^ - pjg-) 

Ap(x) = Ad(x) + Mr) e (\x\ - r) ( - 1 + yfe . 


Lw = d t w — 5Aw 4- b ■ Vw + - (w — (w A d) V (— d)). 

£ 


(36) 


Then, recalling that there exists a positive constant C such that Ad < C and choosing e sufficiently small, 
such that, 1 — yfe > 1 — yfe > \ we have 


Lp- f = -SAd+ M5r) e (\x\ - r)(l - + b- (v d + M yfe rj e (\x\ - r) jfy) + M6 (l - i] e {\x\ - r)) - f 

> - 6 C + M%r) e (\x\ - R) - |b| - \b\My/er] e (\x\ - R) + MS( 1 - rj e (\x\ - r)) - ||/||l~(o,t) 

= M( 6 +(l-\b\y/e- 6 ) Ve (\x\ -r))-SC-\b\ - ||/|Uoo (0 , r) . (37) 

Observe now that the term | — |b| yfe — 5 is negative and, since ?j e (|a:| — r) < 1, we have the following 
inequality 

M (<5 + (f - |b| yfe- S)r] e (\x\ - r)) > M(f - |b| yfe). 

We can fix e 0 such that, for 0 < e < £o, we have |b| yfe < From ( [37] ) , we obtain then 

Ltp - f > M f - S C - |b| - ||/||l~ (0 ,t) = 0, 


provided 

M=^, C' 1 =4(dC+|b| + ||/|| i » (0jT) ), 


(38) 
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concluding then that ip is a supersolution of (34). Analogously, ip is a subsolution of (34) and so we have 
ip < z eS <Tp in Qt and z sS (xo,t) = Tp[x o) = <p[x o). 


(39) 


Observe that, from (36), we obtain 

\Vz eS (x 0 ,t)\ < max{|V^(a;o)|, |V<p(a7 0 )|} < 1 + ^ 

for an arbitrary point xq £ did at any t £ (0, T). We wish to prove that this estimate is true a.e. in 
Qt- Differentiate the first equation of (34) with respect to amultiply it by z% s and sum over k. Setting 
v = \Vz eS \ 2 and noticing that z% 6 k Az% s k = \Av - {z x s kXk ) 2 we get 

-I P 1 P 

-d t v - -Av + -b ■ Vv + - (v - VS £<5 ■ Vz eS ) < 0, 

22 2 — 

being z e& = z eS — ( z eS Ad) V (— d). Using the Cauchy-Schwartz inequality, we obtain 

2 r 

d t v — SAv + b ■ Vr> H - (v— | VS e5 | ) < o. 

£ 

Multiplying the above inequality by (v — (1 + M^/e) 2 ) + and integrating over Q t , we have 

II \(v(t)-(l+M^) 2 ) + \ 2 + 6 [ |V(t> — (1 + My/e) 2 ) + \ 2 
1 Jn JQ t 

+ I b ■ V(u — (1 + My/e) 2 ) + (v — (1 + My/e) 2 ) + + — f (y — \Vz eS \v^)(v — (1 + Mt/e) 2 ) + < 0. (40) 

J Qt £ J Qt 


’Qt 
Since 


/ b • V(v - (1 + M^) 2 )+ (v - (1 + MV~e?) + = 0 


and 


[ (v- |V* e *|«*)(t>- (l + MVi) 2 )+ 

JQt 

= / (y — v^(v — (1 + My/e) 
J tz eS >d\ 


f (y — v^)(v — (1 + My/s) 2 ) + > 0, 

J {z eS <-d} 


from (40) we conclude that (v — (1 + My/e) 2 ) + = 0. 
Then, recalling the choice of M done in (|38|, we have 


|Vz £ | 2 = v < 1 + a.e. in Qt, 

and {z eS } e is uniformly bounded in L°° (0, T; Wq ’°°(fl)). Using (39), it is easy see that 

5. 


—Ci < ~(z eS - (z sS Ad) V (-d)) < Ci. 
£ 


In fact, in the set {z e ° > d} we have 

-(z eS - ( z eS A d) V (-d)) = -(z £l5 - d) < -(Tp- d) < Ci, 

£ £ £ 

in the set {— d < z eS < d} we have z eS — ( z eS Ad) V (—d) = 0 and in the set {z sS < —d} we have 

6 {z eS - (z sS A d) V (-d)) = ^(z £,s + d) > ^ (p + d) > -Ci. 

£ £ £ ~ 


(41) 
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Solutions for linear conservation laws with gradient constraint 


Multiplying the first equation of (34) by d t z eS , we obtain 


f \d t z eS \ 2 + S f Vz eS -Vd t z eS + f b-Vz sS d t z ES + - [ (z sS — (z e6 Ad) V (—d))d t z eS = ( fd t z eS 
JQ.t Jo. t JQt £ JQ.t J Qt 


>Q 
and so 


' Qt 


\d t z‘ 


eS 1 2 


t 2 y n l v ^ tf WI 2 

<U |V^| 2 +(|6|||V^|| i2(QT) + ||J(^-(^Ad)vM)|U 00( Q T) 

1 Jn £ 

+ Wf\\L 2 (0,T))\\dtZ ES \\L2(Q T ) 

|Vu §| 2 + (|b|(l + ^ \fe + Ci)\Q t \ 2 + \\f\\L 2 (o,T))\\dtz E5 \\L 2 (Q T ), 

l^ u o | 2 + 2 (I^IC 1 + T a/e + Ci)\Qt\ 2 + ||/||l 2 ( 0 ,T)) + 2 ll^^lli^Qr) 

where \Qt\ denotes the Lebesgue measure of Qt- So, for <5 fixed, 

\\d t z ES \\ 2 LHQT) <6 [ |V«§| 2 + (|6|(l + ^Vi + C 1 )|Q T |5 + ||/|| L3( o iT )) 2 . 

J Cl 

Then, there exists z s G L°°(0,T; W 0 1 ,oo (O)) n H 1 (0, T; L 2 (fl)) such that 

~ eS v z 6 in L°°( 0 , T; W 0 1 ’ oc (f2))-weak * and d t z eS - A d t z s in L 2 (Qt)- 


< 6 
< 6 


(42) 


£—>0 


£—^0 


Multiplying the first equation of the problem ( |34| ) by v — z sS (t ), where v G K() and integrating over x (s, t), 
0 < s < t < T, we obtain 

d t z eS (v- z eS ) + S I I \7z Ed ■ V(v - z Ed ) + I I b-Vz Ed {v-z Ed ) 

s J Cl 


s J Cl 


s J Cl 


+ - [ (z eS - (z eS (t) Ad) W(-d)){v-z sS ) = N f s {v~z eS ). 

6 J o Js Jn 


For v G Ky, the operator P e defined in (35) is monotone, we have 


[ [ (z eS - {z eS Ad) V {-d)) (■v - z eS ) < 0. 
J s J Cl 


So, letting e —> 0, we obtain 


dtz°(v — z°) + 6 / / V; • (v — z°) + 

J s J Cl 


s J Cl 


b ■ Vz 5 (v - z s ) > 


s J Cl 


f\v-z s ). (43) 


By (39), the function z s is such that z s (t) G Ky, for a.e. t G (0, T). 

To prove that {d t z s }& is bounded in L 2 (Q T ), let e —> 0 in (42), obtaining 


[ \dtz s \ 2 <5 I |Vu 0 | 2 + (|6|(1 + C\)\Qt \ 2 + ||/||l 2 (o,t)) • 

/Qt JO. 


Analogously, letting e —> 0 in (411, we obtain 

|Vz 5 | < 1 a.e. in Qt- 


We can now pass easily to the limit when <5 —► 0 in inequality (43). Observing that z° converges to 
some function z weakly* in L°° (0, T; Wq ’°°(f2)) and d t z s converges weakly in L 2 (Qt) to d t z, we find for all 
0 < s < t < T 


d t z(v — z) + f j b-\7z(v — z)> f f f(v — z) 

J s Cl t/ s */ Cl 
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and it follows also 


f d t z(t)(v — z(t)) + I b ■ V z(t) (v — z(t)) > f f(t)(v — z(t)) for a.e. t G (0, T). 

Jq Jq Jq 

Since z s (t ) G K() for a.e. t G (0 ,T), we also have z{t) G K() and the proof of existence of solution for the 


variational inequality (331 is complete. The uniqueness is also clear. 


The inclusion Ki C K() and the fact z(t) G Ki for a.e. t G (0,T) implies that the function 2 also solves 

□ 


the problem (101. 


Remark 4.2 The first order variational inequalities of obstacle type have been introduced by Bensoussan 
and Lions in m and have been studied in m and in m for general linear operators and general obstacles, 
and extended to a quasilinear two obstacles problem in im ■ In all those cases the notion of solution is less 
regular and the boundary data can only be prescribed on part of the boundary. In addition, the solution cannot 
have a gradient in L 2 and the best that can be expected in general is the operator dtu + b ■ Vit + cu G L 2 , 
as a consequence of Lewy-Stampacchia inequalities. These estimates can be obtained from the regularized 


parabolic inequality (43) and, as in m, it allows the passage to the limit S —> 0 without the estimates on 
the gradient and on the time derivative. It is an open question to establish the equivalence of the first order 
obstacle problem with the variational inequality with gradient constraint for more general first order linear 
operators. 


(44) 

(45) 

(46) 


Theorem 4.3 In addition to the assumptions (|32j), suppose 

b • V^o < f(t) in {x £ Cl : —d(x) < £ 0 ( 2 )} for t > 0 , 

/ = f(t) is increasing and nonnegative, 
liminf/(t) > | 6 | + 2D, 


where D = ||d||£oo(Q) = maxd(r, dLl). Then there exists T* < 00 such that the solution z of the variational 

x(z£l 


inequality (10), or equivalently of (33 1 , satisfies 

z(t) = d for all t > T*. 


Proof We consider z as the solution of the variational inequality (331 . 

Step 1: Zq < z{t) for all t > 0. 

Let v(t) = z{t) + (20 — z {t)) + an d n °t e that v(t) G K(). Then 

[ d t z(t)(z 0 - u(t)) + + I b- L7z(t)(z 0 - z(f)) + > [ f(t){z 0 - z(t)) + . 
Jq Jq Jq, 

On the other hand, by ( |44| , we have 

b ■ Vz 0 (z 0 - z(t)) + < f(t)(z 0 - z(t)) + in {-d<z 0 }, 
and also on {— d = zq} since, in this last set, (20 — z(t)) + = 0 (recall that / > 0). Then 

[ d t z 0 {z 0 - z{t)) + + I b- z 0 (z 0 - z{t))+ < [ f(t.){z 0 - z{f))+. 
Jq Jq Jq 


(47) 


(48) 


From (47) and (481 we get 


f d t {ii 0 - z(t))(z 0 - z{t)) + + / b ■ V(z 0 - z(t))(z 0 - z(t)) + < 0. 

JQ. JQ 


I b-\7(z 0 -z{t))(z 0 -z{t)) + = ]- I b- V((z 0 - z{t)) + f = f \7 ■ b ({z 0 ~ z{t )) + ) 2 = 0 
Jq z Jq z Jq 


But 
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and so 


\ [ | (so -z{t)) + \ 2 <\ [ \(z o -zm + \ 2 =0, 

which implies that zq < z(t), for all t 0 . 

Step 2: z(t) < z(t + h) for all t,h> 0. 


Observe that v(t) = z(t + h) — (z(t) — z(t + h)) G K(), so we can choose v(t) as test function in (331. 
Noting that 

v(t) — z(t) = z(t + h) — z(t ) — (z(t) — z(t + /i)) = ~(z{t) — z(t + /i)) + 


we get 


- f d t z(t)(z(t) - z(t + h)) + - f b-Vz(t)(z(t) - z(t + h)) + > - f f(t)(z(t)-z(t + h)) + . (49) 

Choosing v(t) = z(t + h) + ( z(t ) — z(t + h)) + as test function in ([33]) in the instant t + h and observing 
v(t) — z(t) = z(t + h) — z(t ) + (z(t) — z(t + /i)) + = (z(t) — z(t + h)) , 


that 

we have 


j d t z(t + h)(z(t) — z(t + h)) + f b ■ V z(t + h)(z(t) — z(t + h)) > J f(t + h)(z{t)—z(t + h)) . (50) 


From (49) and (501 we get 


f d t (z(t) - z(t + h))(z{t) - z{t + h)) + I b ■ V (z(t) — z(t + h)) (z(t) - z(t + h)) 

Jn J n 

< [ {f(t)~f(t + h))(z(t)-z{t + h))~ <0, 
J n 

because f(t) < f(t + h), by assumption (45) and (z(t) — z(t + h)) > 0. As 

[ b ■ V (z(t) — z(t + h)) (z(t) — z(t + h)) = 0 , 

J n 


we obtain 


(z(t) - z{t + h)) 


j n (l(^(°) _ z ( h )) 


< o, 


using Step 1. So z(t) < z(t + h), for all t,h> 0. 


Step 3: There exists G %?(S1) such that lim z(x,t ) = Zoo(x), uniformly in x G SI. 

t —>-+oo 

Since the sequence of continuous functions {z(t )} i>0 is increasing in t and is bounded from above by d, 
this conclusion follows immediatly. 

However, in this special case we have a finite time stabilization. 

First we prove that the function Zoo coincides with d. We recall that dtZ G L 2 (Qt), for any T > 0, and 

we set = / z(t). Observe that ||'0||l“(o,oo) < |S1|14, where |fl| denotes the Lebesgue measure of SI. 

Jn 

Since {z(f )} i>0 is increasing, then dt.z > 0 and 


i’(t) ~> Zoo, 

t- S-+00 J n 


> 0 for a.e. t > 0 . 


liminf dt.z(t) =0 in L^fl). 

t—¥ OO 


This implies that 
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Choosing v = d as test function in (33) we obtain, for a.e. t £ (0, oo), 

[ d t z(t)(d - z(t)) + [ b ■ V z(t)(d - z(t)) > f f(t)(d - z(t)) 

JQ JQ JQ 


and so 


d t z(t)(d-z(t)) + \b\ [ (d- z(t))> f(t) f ( d-z(t )). 
! JQ JQ 


Since d > z(t), taking liminf to both sides of the inequality and using the assumption (46), we obtain 

£—>•00 - 


\b\ f (d - Zoo) > (\b\ + 2D) [ ( d-Zoo ), 

JQ JQ 


which is a contradiction unless Zoo = d. 

Consider the following subsets of Q ^ = SI x (0, oo) 


A = {—d<z<d}, I + = {z = d}, I = (z = — d). 


Since z solves the two obstacle problem (|33j), it is well known that the following inequalities are verified 
a.e. in 

dtz + b-\7z = f in A, dtz + b ■ V z < f in I + , dtU + b-Vu>f in I~. 

If there is no finite time stabilization of the solution, since z(t) is increasing in time, we may find a point 
(xo, to) an d an open subset oj o of S7 with xo £ tc 0 , such that, (x, t) £ A U I~ for t > t 0 . So, 

/(f) < d t z[x , t) + b ■ Vz(x, t) for a.e. (x, t) £ wo x [f 0 , +oo). 

Then, for any t > to and any open set u C wo, we have 

rt+1 pt+1 /* 

f{r)<j^Jt j (dtz(x,T) + b-Vz(x,T)) 

< 1^1 / (z(x, t + 1) - z(x, t)) + \b\ < (2D + |6|). 


As a consequence, 


rt+1 

liminf/(t) < liminf / f(r)dr < 2D + \b\ 

£—>■00 £—>• oo J f 


and this is a contradiction with (46). So z(t) must stabilize in finite time. 


□ 
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